When hot water freezes before cold 
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^ ■ Abstract 

^\j ■ I suggest that the origin of the Mpemba effect (the freezing of hot water before cold) is freezing- 



point depression by solutes, either gaseous or solid, whose solubility decreases with increasing 
temperature so that they are removed when water is heated. They are concentrated ahead of the 
-.^ freezing front by zone refining in water that has not been heated, reduce the temperature of the 

Q ■ freezing front, and thereby reduce the temperature gradient and heat flux, slowing the progress of 

Q \ the freezing front. I present a simple calculation of this effect, and suggest experiments to test this 

^' hypothesis. 
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In a sub-freezing environment initially hot water often freezes before initially cold water. 
This observation is counter-intuitive because one naively expects the hot water first to cool 
to the temperature of the initially cold water, and then to follow the cooling history of the 
initially cold water. However, the effect has been observed many times and is folk-wisdom in 
many cultures; the earliest known reference is by Aristotle. It was brought to the attention 
of modern science by the Tanzanian high school student for whom it is now named; with 
admirable persistence in the face of disbelief on the part of his teachers, he insisted on the 
primacy of empirical evidence over theory. The history and literature are summarized by 
Auerbach |l|. 

No generally accepted explanation of the Mpemba effect exists. Apparently, pre-heating 
water affects its properties in a manner that accelerates its freezing. A number of mechanisms 
have been considered, including loss of mass by evaporation '^, the loss of dissolved gases 
(whose solubility in hot water is much less than in cold water) and supercooling 1]. 

Wojciechowski, .i at. B .eport what appear to be the only systematic cuantitative 
measurements of the Mpemba effect. Following a suggestion of Freeman |^, they also 
measured the freezing of water saturated with CO2. Unfortunately, they did not describe 
the ionic content of their water, which is likely to be essential. 

Auerbach |l| found substantial (several degrees) but non-reproducible supercooling in 
both preheated and non-preheated samples. As a result, some of his data showed an Mpemba 
effect while some did not. His experiments were performed on distilled and de-gassed water, 
and may therefore not be applicable to observations of the Mpemba effect in tap water or 
environmental water (Mpemba's observations were on sugared milk he was freezing to make 
ice cream!). 

The observations clearly point to some change in water when heated. As has been re- 
marked before, heating water removes dissolved air (chiefly nitrogen) because its solubility 
decreases rapidly with increasing temperature. The problem is to find a mechanism by which 
the removal of a small quantity of dissolved material (the solubility of nitrogen in water at 
room temperature is only about 0.7 mmolar/bar) can produce a Mpemba effect. 

Gases are not the only substances whose solubility in water decreases with increasing 
temperature. Most natural waters are "hard", containing a variety of dissolved mineral 
salts, most importantly calcium bicarbonate Ca(HC03)2. This is introduced into ground 
water that has been acidified by atmospheric carbon dioxide by the reaction in limestone 
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FIG. 1: Temperature profiles, neglecting specific heat in comparison to latent heat. Tm is the 
freezing point of pure water, T^ < T^ that of hard water, lowered by freezing point depression. Its 
temperature gradient and heat flux are reduced, so the freezing front advances more slowly than 
in water purified by heating. 

rock [5| 

CaCOs + CO2 + H2O ^ Ca(HC03)2. (1) 

Because the solubility of gases in liquids decreases rapidly with increasing temperature, 
heating shifts the equilibrium to the left, resulting in the precipitation of limestone deposits 
known as "boiler scale" or "kettle fur" . For this reason, hardness resulting from bicarbonates 
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dissolved calcium. Hence the freezing point T^ of never-heated hard water is lower, because 
of the depression of freezing points by solutes |6] , than the freezing point T^ of hard water 
that has been heated. For the same thickness of the ice layer, the temperature gradient 
and heat flux are less in the never-heated water, so it takes longer to lose the latent heat of 
freezing and freezes more slowly. This is illustrated in Fig. 1. 

The concentration of dissolved calcium in moderately hard tap water is about 100 ppm, or 
2.5 mmolar, but values several times greater or less are found. Each Ca"'""'' ion corresponds to 
a total of three dissolved ions, so such a solution has a total ion concentration of 7.5 mmolar. 



The freezing point depression of ideal dilute aqueous solutions [6|] is ATm = RT^/Hm = 
1.86° C/molar, where R is the gas constant and Hm the enthalpy of melting. The effect 
of hardness is typically ten times greater than that of dissolved air, but still rather small 
without further concentration. 



a solute in a solid to that in its melt is the 
. For most solutes in water and ice fco ^ 1; 



The equilibrium ratio of the concentration o^ 
Nernst equilibrium distribution coefficient fco 
they are almost completely excluded from the solid. As a freezing front advances into water 
the solutes are pushed ahead of it, and their concentration close to the front is enhanced 
over its initial value in the liquid, a process known as zone refining |2|. 

The conventional theory of zone refining |S] assumes a stationaryi state with a boundary 
condition that the concentration of solute equals its initial concentration in the liquid at 
a finite distance from the freezing front. This latter condition is appropriate when there 
is driven fiuid circulation outside a viscous boundary layer of thickness 6. In the present 
problem there is no circulation and 6 -^ oo. Worse, in the limits ko ^ and 6 ^ oo this 
solution becomes singular because in steady state the solute from an infinite volume of fiuid 
must accumulate in a finite region of enhanced concentration. 

Our problem corresponds to the low concentration limit of the theory of the freezing 
of binary melts '0\, in which the progress of the freezing front is controlled by diffusion of 
heat. Mathematically, it is a "Stefan problem" in which the boundary condition is moving. 
We neglect the specific heats of the liquid and solid in comparison to the latent heat of 
the phase transition. Then the liquid is isothermal and there is no thermal diffusion of the 
solute within it. 

The freezing front advances at a nonsteady speed f/(t) = J{Tjn — Tq)k/ {2HmPit) |lO|, 
where Tm is the melting temperature, Tq < Tm is the subfreezing boundary temperature, k 
is the thermal conductivity and pi the density of ice, and t is the time elapsed since freezing 
began. The equation for the conservation of solute concentration C{x,t) in the liquid, in a 
frame moving with the freezing front, is 

dt +''^*^ dx dx^ "^' ^^ 

where v{t) = —Vf{t)pi/pi is the velocity at which liquid flows to the freezing front and D is 
the diffusion coefficient of solute in the liquid. 

Both the diffusion length y/Dt and the frozen layer thickness J Vf{t) dt are proportional 



to t^/^. Hence the solution is stationary in suitable variables. Defining ^ = x/yDt juj, we 
obtain the separable equation 



where the dimensionless parameter Km = J{Tm — To)Kpi/{2HmPiD). The value of D 
for Ca(HC03)2 in water is not readily available, so we take the values |lj, |l3| for 
CaCl2 as estimates and extrapolate to 0°C, giving D ^ 0.6 x 10^^ cm^/sec. Then 



i^M ~7^(r^-To)/10°C. 

Writing C{(,t) = C^(^)Ct(t), we find ordinary differential equations for C^(C) and Ct{t) 
with a separation constant A: 

dlnCtit) 
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A (4) 



-^^ + {Km + -)^^ - ACdO = 0. (5) 

For Km ^ 1 and ( = 0{1) we can neglect the (/2 term and the solution is elementary. From 
the conservation of solute (and ko — * 0) we find A = 0. The equation for C(^(C) has solutions 
C(^{() oc exp (7^) with 7 = 0, —Km- The root 7 = gives the uniform solute density far 
from the freezing front, while the root 7 = —Km gives the enhanced concentration of solute 
close to the front: 

Cc(C) = Co[l + 2Kl,^exp (-KmC)]- (6) 

Pi 



This justifies the neglect of the (/2 term to 0{K 
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The concentration at the freezing front is enhanced by a factor ^ 100(Tm — To)/10°C. 
The fractional freezing point depression 

where Cq is the concentration in moles/cm^, independent of time and of Tq. For Ca(HC03)2 
in water the coefficient of Cq is 20/molar so an ion concentration of 5 mmolar (67 ppm 
calcium) produces a 10% reduction in freezing rates compared to pure (once heated) water. 
Wojciechowski, et al. |j| reported an apparent increase of Hm of C02-saturated water of 
10 ± 4 cal/gm. This is not likely to be an actual change in Hm-, but may refiect a reduction 
of the freezing point at the freezing front of 10 ± 4°C, requiring the removal of additional 
internal energy before freezing begins. Application of the theory presented here predicts 



that for water saturated with respect to CO2 at 1 bar pressure at room temperature (33 



mmolar CO2) and Tm — Tq = 23 °C 



a 



the freezing point depression at the front is 10 °C, 
consistent with the apparent increase in Hm- These authors found an Mpemba effect only 
in water that had not been degassed by boihng, which is consistent with the removal of 
temporary hardness, even though dissolved air itself is predicted to produce only a small 
Mpemba effect. 

Our model makes readily testable predictions for the dependence of a Mpemba effect 
on the concentration of solutes whose solubility decreases with increasing temperature. It 
also predicts no Mpemba effect, but a dependence of freezing time on concentration, for 
solutes whose solubility does not decrease with increasing temperature. It predicts that the 
magnitude of freezing point depression and of the effect (having subtracted the time required 
for the solutions to cool to the freezing point) should be independent of the thickness of ice 
formed. In this model an Mpemba effect occurs in the time required for formation of an ice 
layer of any thickness, and hence is found for the onset of freezing as well as for complete 
solidification. Fully quantitative predictions require modeling of heat transfer through the 
air boundary layer. 
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